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ABSTRACT 

Measurements of Pu-Be neutron age to indium resonance 
and thermal migration area were made in pure water at 2 ^°C 
and in mixtures containing 1 . 16 , 8 . 96 , 27 . 81 , 46 . 65 3 
70 . 09 , and 80.09 grams/liter of boric acid at 40°C. Thermal 
diffusion lengths were calculated from these measurements 
and the dependence on poison concentration observed. It 
was concluded, both from failure of the age measurement 
and the irresolute thermal flux behavior at the heaviest 
poison concentration, that the "maximum absorption" limit 
was exceeded in this experiment. The absorption cross- 
section of the supporting medium corresponding to the 80.09 
g/1 concentration was 0.6120 cm“^. The present results, 
together with previous experimental data, provide an esti- 
mate of diffusion length at the absorption limit of 0 . 5 cm, 
slightly less than that corresponding to the theoretical 
Corngold limit (0.6 cm). 
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I . INTRODUCTION 



It has been shown analytically by Corngold [1] that if 
the absorption concentration of a supporting medium is 
increased beyond a certain value (assuming 1/v absorption 
in an infinite medium), no thermal neutron diffusion length 
exists. The practical implication of this is that under 
such extreme circumstances the spatial distribution of 
thermal neutrons cannot be resolved. Table 1 lists theo- 
retical values of the minimum diffusion length which 
corresponds to the "maximum absorption" limit for the more 
common moderators [2] . 



TABLE 1 


Moderator 


Minimum Diffusion 
Length (cm) 


C 


23 


Be 


18 


o 

C\J 


0.6 


D^O 


2.8 



The dependence of thermal neutron diffusion length upon 
the absorption concentration in water has been the subject 
of several experiments [3> 53. In all of these experi- 

ments the limiting value of diffusion length is assumed, 
but not proved. 
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The most extensive measurements of diffusion length In 
water poisoned with a 1/v absorber have been those by 
Goddard and Johnson [33. Their work Included measurements 
In water containing dissolved boron which was Ideally 
suited for the purpose. Boron Is a strong 1/v absorber, 
l.e.jlts cross-section varies Inversely with neutron velo- 
city, or as the square root of neutron kinetic energy. 
Therefore, the homogenlous mixing of small amounts with 
pure water has negligible effect on the average scattering 
cross-section but drastically Increases the absorption 
cross-section . 

Goddard and Johnson experimented with a wide range of 
boron concentrations at temperatures of 22.3° and 65°C. 

Their most heavily absorbing solution contained 69*3 g/1 
of boric acid (H^BO^) and gave a diffusion length measure- 
ment of 0.6581 cm (within the Corngold limit). Similar 
measurements using dissolved cadmium at concentrations of 
55.5 and 82.2 g/1 resulted in measurements of 0.5983 and 
0,4925 cm respectively, both of which lay outside the 
Corngold limit. They did not report any difficulties In 
obtaining these latter diffusion lengths. 

This research project aimed at experimentally determin- 
ing the "maximum absorption" limit by Increasing the boron 
concentration In pure water beyond that attained by Goddard 
and Johnson until the diffusion length measurement failed. 

This experiment was conducted using a plutonium— beryllium 
source which emitted 5.^ x 10^ fast neutrons per second. 
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It will be explained in the following sections why this 
fast neutron source dictated an unconventional and Indirect 
method of determining diffusion length. In the process, 
information concerning the dependence of neutron migration 
area and age on the poison concentration is obtained for 
the first time. 
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II. THEORY 



A. TWO-GROUP DIFFUSION THEORY 

The thermal neutron flux, in most cases, arises from a 
source of neutrons having energies many orders of magnitude 
larger than the energy of neutrons in thermal equilibrium 
with the medium (about .025 eV at 25°C) [6]. The Pu-Be 
(a,n) source used in this work yields neutrons with an 
average kinetic energy of about k.2 MeV and the neutron 
spectrum ranges from roughly zero to about 10.6 MeV. The 
calculated spectrum is shown in Figure 1. (7l 




The analytical model which describes the neutron popula- 
tion as two separate groups, one fast and one thermal, gives 
solutions which approximate the behavior of neutrons in real 
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assemblies and for many practical calculations. Is suffi- 
ciently accurate [8]. In the two-group treatment, neutrons 
of thermal energy constitute one group vfhlle all those of 
higher energy belong to the fast group. Neutrons are 
considered removed from the fast group by scattering 
collisions which reduce their energy to the thermal region. 
Since absorption cross-sections of moderators for high 
energy neutrons are very small, any loss of the fast 
neutrons by absorption Is neglected. 

The disappearance of neutrons from the fast group Is 
described In terms of a "slowing down" cross section, Ip 

16] . By denoting the fast neutron flux by <j)p (neutrons/ 

2 

cm -sec). Ip Is defined such that Ip4)p(r) gives the number 

o 

of neutrons per cm /sec. which slow down out of the fast 
group at the point r cm from the source. 

If an isotropic point source at r = 0 emits S fast 
neutrons/sec in an infinite moderator, the fast flux at 
r / 0 can be described by [73 

DpV^(J)p(r) - Ip4>p(r) = 0 (1) 

2 

where V 4>(r) is the Laplaclan of the flux and Dp is the 
diffusion coefficient for the fast neutrons. The source 
at r = 0 is taken into account by means of a source con- 
dition which will be stated later in this section. 

The quantity Dp can be derived by considering Pick's 
Law of Diffusion which states that the neutron current 
density Is proportional to the negative gradient of the 
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flux. The proportionality constant has dimensions in 
length and is called the diffusion coefficient, denoted 
by D in the more general treatment of neutron diffusion. 
Accordingly, the fast neutron current density can be 
represented by [6] . 



Fast Neutron 
Current Density 



= - DpV(|)p(r) = - / D(E)V(j)(r,E) dE 

^th 



where the Integral is evaluated over the energy region of 
the fast froup, i.e., from thermal, to source energy, 

Eq. The fast flux is defined by the Integral 



E 

(t>p(r) = / ^ i()(r,E) dE 



and therefore 

^0 

/ D(E) V({)(r ,E)dE 
E4.U 

D = -th 

E 

/ ° V(j)(r,E) dE 
Eth 

In two-group theory it is assumed that the variation in the 
neutron flux spectrum from point to point in the medium can 
be neglected, that is the flux is completely separable in 
space and energy. Thus, 
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D„ = 



0 

"th 



/ " D(E)(f)(E)dE 

E. 



( 2 ) 



/ <1)(E) dE 

^th 



where <})(E) is the energy dependent flux. 

Returning to Equation (l), and using the boundary con- 
dition that <f)p(r) must remain finite as r-»-“, the general 
solution can be obtained: 



-r//r^ 



<}>p(r) = C - 



where the parameter Is known as the neutron age to 
thermal energies, defined here as [9] 



( 3 ) 



D. 

^ Z, 



P 



(^) 



The constant C Is found from the source condition which is 

obtained by drawing a small sphere around the source and 

computing the number of neutrons that pass through Its 

surface per second. Considering the magnitude of the cur- 

2 

rent density In fast neutrons per cm /sec, Jp(r) , and a 

2 

sphere of radius r, this number Is ^Ttr Jp(r). Then, in 
the limit, as r 0 , the source condition Is 



lim ^lirr J„(r) = S 
r -^0 ^ 



( 5 ) 



13 



Thus 



-rZ/xT 



T 



d<l)p(r) 





r 



2 



and 




F 



Finally, the fast flux is given by 




( 6 ) 



Considering now the thermal neutron group, the time- 
independent diffusion equation which describes the thermal 
flux can be written 



where <I),p(r) is the thermal , flux, D is the average value of 
the diffusion coefficient in the thermal region, is 
the average macroscopic thermal absorption cross section 
and q,j(r) is the thermal slowing down density. The quantity 
qrp(r) is equivalent to the production rate of thermal 
neutrons per unit volume and, for the case of a point source 
of fast neutrons in an Infinite moderator, is equal to 
Ep<|)p(r), the number of neutrons slowing down out of the fast 
group . 

Equation (7) can be rewritten as 



DV^ <l),p(r) - <|),p(r) + q,p(r) = 0 



(7) 



<f>rp(r) - 







0 



( 8 ) 



where <}>p(r) Is given by Equation (2) and = 



D l"" . 

r' 



the thermal neutron diffusion length. By requiring the 
neutron flux to remain finite, the homogeneous solution 
becomes 

-r/L„ 

e ^ 

4>ip(r)H = - . (9) 



Since the source term is proportional to the fast flux, a 
particular solution is assumed of the form 



4>,p(r)p 



K. 



-r//r^ 

e 



3 

r 



(10) 



By substituting this expression into Equation (8), the 
constant Kp is found to be 



K 



2 



-S 

Hit D (L^ - T^) 



( 11 ) 



The general solution to Equation (8) can now be written 

~r/Lrp — ~r/ v^Tfp 

S e 

( 12 ) 



4>Y(r) = 



Kp e 



ilTTr D (L^ - T^) 



where the constant can be determined by a source con- 
dition as in the fast flux case. However, since the source 
only produces fast neutrons, the required source condition 
for the thermal flux is 

lim Jm(r) = 0 (13) 

r^O 
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where J,p(r) is the magnitude of the thermal neutron current 
density. Thus, 



Jrp(r) = -D 
and the constant 






d<}),p(r) 

dr 

becomes 

2 




(l'^) 



Finally, the thermal flux is given by 






ilTir Z^L^2(1 -t^) 



~r/Lrp r/ v^Tm 
(e -e M (15) 



Equation (15) describes the thermal neutron flux in 
2 

neutrons per cm /sec In an Infinite moderator using two- 
group theory . 

It follows from the development of Equation (15) that 

it is also possible to describe the slowing down of the 

fast neutrons and their subsequent diffusion as thermal 

neutrons as two distinct processes. The neutron age, tip, 

is related to the distance traveled from the point where 

a fast neutron first appears to the point where it slows 

down to thermal energies. Similarly, the square of the 

2 

thermal neutron diffusion length, , corresponds to the 
distance traveled from the point where a neutron becomes 
thermal to the point where it is finally absorbed. A 
quantity M,p , known as the thermal migration area, takes 
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into account the total neutron travel during moderation 
and diffusion. Accordingly, 

( 16 ) 

where (migration length), /x^ (slowing down length) and 
L,p all have dimensions in length. 

B. ONE-GROUP DIFFUSION APPROXIMATION 

Although the two-group theory gives more accurate flux 
solutions, the one-group approximation readily accommodates 
the "first scattering source" theory which is useful for 
our analysis. In one-group diffusion theory all neutrons 
are lumped together into one-group, namely thermal. This 
approximation, while limiting the description of the neutron 
population to spatial dependence, allows attention to be 
focused on the diffusion process. 

In the two-group model of the previous section, neutrons 
entered the thermal group as a result of slowing down out 
of the fast group, and the thermal slowing down density, 
q,j,(r), appeared as the source term in the thermal diffusion 
equation. In the one-group treatment, however, when the 
actual neutron source yields fast neutrons, the source term 
used for the diffusion equation must be altered to account 
for the behavior of the fast neutrons prior to reaching 
thermal energies. In the particular case of a point source 
emitting S fast neutrons/sec in an infinite media, the 
slowing down density at r is given by the relation [10] 
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q(r,T) = 



-r^/^T 



(17) 



Se 

where t Is the neutron age. A small t means that the neu- 
trons have suffered little slowing down and so they have 
. not diffused far from the source. Larger values of t 

correspond to neutrons which have been slowing down longer 
and therefore have had an opportunity to diffuse further 
from the source. When the neutron age to thermal energies, 
Ty, is substituted into Equation (17), the resulting expres- 
sion gives the dispersion of thermal neutrons and can then 
be used as the source term for the diffusion equation. 

Equation (17) is based on the Fermi age model which 
assumes that the slowing down process is a continuous one 
[9). In hydrogeneous media, however, the slowing down 
process is not continuous due to the fact that a neutron 
may lose all of its energy following a collision with 
hydrogen. In addition, the hydrogen scattering cross 
section decreases rapidly with energy above about 10 keV. 

As a result, the fast neutrons found far from the source 
are those resulting from high energy neutrons which have 
not suffered a collision. When these neutrons finally do 
collide they will be moderated without moving very far from 
the point of collision. Accordingly, a fair representation 

of the source of thermal neutrons which can be used as the 

/ 

source term in the diffusion equation is the "first 
scattering source" [10] which for a point source is given by 
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I 

I 




This equation gives the density of first scattering colli- 
sions at distance r from an isotropic point source emitting 
S fast neutrons/sec in an infinite medium. The parameter 
is the macroscopic scattering cross section of the 
fast source neutrons. The reciprocal of is known as 
the relaxation length, b, defined as the distance the 
neutron density is reduced by a factor 1/e of its Initial 
value due to scattering. Replacing E^ by 1/b , the source 
term can be written 



At large r, therefore, the principal variation in the 
slowing down density arises from a simple exponential 
decrease combined with the inverse square characteristic 
of a point source. Moreover, q(r) is not a function of 
3-Se > since every neutron is presumed to slow down 

fully at its first collision. 

The thermal diffusion equation for the one-group approx- 
imation can now be written as 



q(r) 




(18) 



5 

where “ ipff ^ constant 




-r/b 



2 



0 



(19) 



Dr 
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To satisfy the boundary condition that 4>,p(r) go to zero as 



r -»• oo, the homogeneous solution becomes 



-r/Ln 






Ke 



( 20 ) 



A particular solution is then assumed of the form 

° — ~:2 • 



( 21 ) 



By substituting this expression into Equation (19) > the 
constants and are found to be 



= 0 



and 



^2 = 



2 2 
2 2 ” 



( 22 ) 



The general solution to Equation (19) can now be written 

(23) 



<}>,p(r) = 



K e 






Dr^(L^^ - b^) 



If the diffusion length, L^, is larger than the relaxa- 
tion length, b, the first term of Equation (23) dominates 
at sufficiently large r, the thermal flux decays expo- 
nentially, and can be evaluated in a simple routine 
manner. This is the conventional method of determining Lr^. 

The relaxation length, however, depends on the initial 
energies of the neutrons generated by the source. Fast 
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i 



neutrons from a Pu-Be source have relaxation lengths of 
about 10 cm. in water [11] which exceed 2.78 cm., the cor- 
responding diffusion length of thermal neutrons [6] . Con- 
squently, the thermal flux decays In a more complicated way 
and the simple theory for the case of a completely thermal- 
Ized neutron field does not apply. The determination of 
In water can, however, still be accomplished by measuring 
both the migration area and the neutron age. From Equation 
(l6) the diffusion length is then obtained from the relation 

Ly = 

C. DETERMINATION OP NEUTRON AGE 

The neutron age, x, between two energies Is proportional 

to the mean square distance travelled by a neutron as it Is 

being slowed down from the higher to the lower energy. 

_ 2 

Denoting this distance by r^ , the second spatial moment of 
the slowing down density can be written as [93 



rg 2 (Eg,E) 



/ r^q(r, T(Ef.,E) )dv 
/ q(r,T(EQ,E))dv " 6 t(Eq,E) 



(25) 



where Eq is the energy of the neutron at Its point of birth 
as a fast neutron and E Is a lower energy attained as a 
result of the slowing do\vn process; dv Is the volume element, 
and q is given by Equation (17). Neutron age Is now defined 
by the relation [93 

F 2 

t(Eq,E) = (Eq,E) (26) 



Equation (25) applies to neutrons of any energy and. 



hence, of any age. In actual measurements, however, the 



terminal age Is fixed by the resonance energy of the detect- 
ing foil. Consequently, in determining the age to thermal 
energies, x,j,, the spatial distribution of the slowing down 
density at some energy above thermal is obtained and then a 
correction is applied to account for the additional age 
required to reach the thermal region. 

The slowing down density of source neutrons at 1.^5 eV 
can be determined by activating cadmium-covered indium foils 

placed at various distances from the source. The principal 

115 

indium isotope, IN , has a strong. Isolated absorption 
resonance at 1.^5 eV while the absorption cross-section of 
cadmium is large below 1.^5 eV, but small at higher energies. 
The cadmium cover, therefore, prevents thermal neutrons 
from reaching the indium foil; whereas the large indium 
resonance at 1.^5 eV is primarily responsible for any 
neutron induced activation. 

The saturation activity of a cadmium covered indium 
foil, is equal to the rate of production of 

nuclei and can be written as [9] 

A_„ = V / (})(r,E) (E) dE (27) 

CD ajj^ 

where V is the volume of the foil, E is the macroscopic 

115 

absorption cross-section of IN and the Integral is 

evaluated across the 1,^5 eV resonance. The number of 
neutrons per cm /sec which slow down past a given energy, 

E, in the diffusion approximation, is given by 

q(r,E) = <l)(r,E)C I^(E)E (28) 

where <})Cr,E) is the energy dependent neutron flux, E (E) 

o 

is the macroscopic scattering cross-section, and 5 is the 



average logarithmic energy decrement. From Equation (28) 
it follows that 



^CD = V 



rq(r,E) Eajjj(E) 



vet 



E 



(29) 



Assuming that q(r,E), and 2^^(E) are constant over the 
small energy range of the resonance peak, this can be 
written 



^ Vq(r,1.45) ^ j 

CD C E ; 

s 



•IN 



(E) 



dE 



(30) 



That is, the slowing down density at 1.^5 eV is proportional 

to the experimentally measured activity of the foil. By 

substituting this result into Equation (25) j and noting 
2 

that dv = ilTrr dr, the expression for the age to indium 
resonance in spherical coordinates is 



^IN 



- 2 



/ r App, dr 
10 



/ r^A^^dr 



CD 



(31) 



Thus, by knowing only the relative values of q(r, 1.^15) j 
the age to 1.^5 eV can be obtained by numerically evaluating 
the integrals in Equation (31). In this analysis the slowing 
down density is assumed to be approximated by Equation (l8). 

Another method for measuring neutron age [6] , which 
provides only a rough estimate in hydrogeneous media, is 
based on the spatial distribution of neutrons of a given 



energy described by Equation (17). Taking the logarithms 
of this equation and considering Indium resonance energy 
leads to 



2 

Hn q(r) = const - r . 



(32) 



Since the resonance activation Is proportional to q(r), a 
plot of £n App measured at various points from the source 
can be made. For moderators which satisfy the continuous 
slowing down theory of the Fermi Age model (e.g. graphite) 
the experimental points will lie on a straight line. For 
water, however, the points will lie above a straight line 
except In a very narrow band about 20 cm. from the source 
[12] . This Is due to the unusual nature of the hydrogen 
scattering cross-section and the fact that Individual 
neutron energy losses are so large, as previously discussed. 
In this narrow band, however, an approximate straight line 
can be drawn whose negative slope will provide an estimate 



D. DETERMINATION OF THERMAL MIGRATION AREA 

The thermal migration area, defined by Equation (l6), 
correlates with the total mean square distance a neutron 
travels from Its point of birth to Its point of capture. 
This Is given by [11] 



of T 



IN 




/ Aq r dr 
0 ^ 



( 33 ) 




2k 



where Ag Is the total saturated activity measured as a 
function of distance r from a point source. Equation ( 33 ) 
can be obtained in a manner analogous to the derivation 
of Equation ( 31 ) for rg . Here, however, it is the total 
relative neutron flux which must be known to evaluate the 
integrals. At energies that activate a bare Indium foil, 
this total flux will be proportional to the saturation 
activity. Migration area can then be determined from the 
definition * 



M 



T 



2 




( 3 ^) 
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III. DESCRIPTION OF THE EXPERIMENT 



Measurements were made In pure H^O at 24 °C and In 
mixtures containing l.l6, 8 . 96 , 27.81, 46.65, 70.09 and 
80.09 grams of H^BO^ per liter of H 2 O at 40°C. The last 
two mixtures contain the highest boron concentrations 
attained to date. 

The experiments were conducted in a cylindrical plexi- 
glass tank, 1.45 ft. in diameter and 4.27 ft. high, which 
held about 50 gallons of water. This tank was situated 
atop a smaller, steel tank 2.95 ft. in diameter and 1.25 
ft. high, which held about 65 gallons of water. The base 
tank served to contain the Pu-Be source and also provide a 
protective water shield. When positioned for measurements 
in the base tank, the source was 0.625 inches below the 
base of the plexiglass tank at centerline. A polyethylene 
Jacket 0.125 in. thick made the source holder watertight. 

A single steel rod immersed in the tank at the side 
heated the poisoned mixtures through a rheostat control. 

This heating rod, plus a continuous running electric stirrer, 
maintained the medium at a uniform temperature of 104°F 
(±0,5°). 

The relative flux was measured with bare Indium foils, 
1.0 in. in diameter and .010 in. thick. Foils used in any 
given measurement had a weight spread of less than 2.65^. 
Indium resonance activity was measured with the same foils 
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placed In cadmium covers 0.025 in. thick. A plexiglass 
foil holder arrangement provided accurate spacing of the 
foils in the medium on a line perpendicular to the source. 
The distance between adjacent foils was 0.50 Inches for 
all measurements. Figure 2 shows the experimental appar- 
atus and set-up. 

Nine foil positions were used for each measurement. 

The closest measuring position to the source was 8.26 cm. 
while that furthest away was 18.^11 cm. Counts of the cad- 
mium covered foils became insignificant beyond this last 
position . 

The cadmium-covered and bare Indium foils were irradi- 
ated separately during each experiment for a minimum of 8 
hours. After irradiation the indium was allowed to decay 
for 3 minutes to avoid interference with the desired 
5^-mlnute activity from the 13 second half life of 
[11] . The activity of the foils was then counted from each 
side for 10 minutes; the relatively long counting time 
having been chosen to Improve statistics. The counting 
system consisted of an end-window Geiger Mueller tube con- 
tained in a lead shielded housing, and an electronic 
scaling unit. This apparatus is summarized in Table 2 
and pictured in Figure 3« 

Counter background was measured repeatedly during each 
run as an erratic line voltage perturbation appeared at 
infrequent Intervals. Only rarely was the background 
stable within statistics for an entire run. The resulting 
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Experimental Apparatus 



Figure 2 
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Counting Equipment 



TABLE 2 



DETECTION EQUIPMENT 



Counter 

Activity Detected 

Type 

Shape 

Counting Gas 
Manufacturer 

Amplifier 

Type 

Sensitivity 

Manufacturer 

Model 

Scaler-Timer 

Manufacturer 

Model 



Beta and Gamma 
Organically-quenched GM 
Cylindrical 
Helium 

Atomic Accessories, Inc. 

Non-overloading 
10 mV 

Hewlett Packard 

555^A 

Hewlett Packard 
5201L 



necessity for interpolating between background determina- 
tions represents one Important possibility for nonrandom 
errors in the results. 

Prior to each poison experiment samples of the solution 
were chemically analyzed for boric acid content. The method 
of analysis was boric acid titration with glycerine which 
is outlined in Appendix A. 
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IV. ANALYSIS OF DATA 



A. DETERMINATION OP SATURATED ACTIVITY 

Three corrections were applied to the foil data. The 
first was an internormalization based on foil weight, the 
second corrected for absorption of indium resonance neutrons 
by the cadmium cover; and the third corrected for counter 
resolving time. 

The weights of the indium foils used in all measurements 
varied from 0.9282 to O .9621 grams. All foil weights were 
therefore normalized to a weight of 1.0 gram and the result- 
ing corrections applied to the foil counts. The highest 
value of the weight correction factors obtained In this 
manner was I.O 82 . 

The Indium resonance counting-rate correction was 
determined to be 1.11, This value was obtained from a 
curve of correction factor vs. indium foil thickness con- 
structed for a cadmium thickness of 0,025 In. as outlined 
by Tittle [1^] . 

A number of determinations of counter resolving times 

116 

were made by following the decay of 5^ minute IN"^ over 
three hours. An average value of 400psec. was obtained by 
solving 



1 





2 

N2t 
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where A Is the extrapolated activity, N Is the mean counting 
rate, and t is the resolving time. A correction factor of 
1/(1-Nt) was then applied to all counting runs. 

The effect of foil activation by source neutrons above 
resonance was considered negligible as the closest measure- 
ment point was greater than 8 cm. from the source [12] . 

Also, to Insure that the close spacing between foils did 
not contribute error due to shadowing, activations were 
made in pure water using various spacing intervals. 

After correcting the foil count as described above, 
the saturated activity In counts per minute was calculated 
from the equation 



A 



S 



( 1 - 




XC 

-X(t2-t^) 



-X(t.-t, ) 
e 3 ^ ) 



where C is the total count less background, X is the disin- 
tegration constant of t^ is the irradiation time, 

^2 ^1 waiting time and t^-t^ is the waiting time 

plus counting time. 1133. An average quantity was then 
obtained from the values of Ag computed for each foil side. 
The computer program used in this calculation is included 
in Appendix P. Appendix I provides tables of the computed 
data for all runs. 



B. DETERMINATION OF RELAXATION LENGTH 

It was pointed out in Section II-B that in hydrogeneous 
media, the slowing down density at large distances from the 
source can be approximated by the expression 
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q(r) 



e 



-r/b 



( 18 ) 



where b Is the relaxation length of the fast source neutrons. 
Since In the framework of age-diffusion theory, the flux 
is proportional to the slowing down density at a fixed 
energy, it follows that the cadmium-covered activity at 
large distances from the source would have the same space 
dependence described by Equation (l8), namely 



Acd(i’) 



K e 
r 



-r/b 

2 



(35) 



To evaluate b from the corrected data, the product 
2 

Acd^’ was plotted against r on semi-log paper. This plot 
showed the spherical space distribution beginning to 
decrease exponentially at about r = 12.5 cm. The distribu- 
tion for r > 12.5 cm was then fitted, by least squares, 
with functions of the form of Equation (35); specifically 

2 

i-n Apj^r = const - r/b (36) 

The distribution for r > 12.5 cm consisted of only five data 
points which, for most runs, exhibited large statistical 
variation due primarily to the low foil activity in this 
range. Accordingly, the solution resulting from fitting 
( 36 ) to this data provided only a rough value of the param- 
eter b. A more accurate determination was made by comparing 
A„„ distributions calculated from Equation (35) (using 
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neighboring values of b) with the measured distribution on 

a large scale semi-log plot of vs. r. Unlike the 

2 

A^^r spherical space distribution all data on the A^^^ vs. 
r plot decreased in value with increasing r; this permitted 
use of the more reliable data close to the source to indi- 
cate which, if any, of the data at r > 12.5 cm were incon- 
sistent with the decreasing trend of the overall distribu- 
tion. After discarding any poor data Indicated by this 
procedure, the calculated A^^^ distributions were normalized 
and fitted to the measured distribution for r > 12.5 cm 
until the closest possible fit was attained. The value of 
b used in calculating the distribution found to give the 

closest fit became, then, the negative reciprocal of the 

2 

slope of the exponential portion of theA^^r vs. r curve. 

The above procedure was also used to determine relax- 
ation length for the bare foil activity (Ag) and, thus, the 

2 

slope of the A^r vs. r curve. Foster, in measuring the 
age of Na-Be neutrons [15] , assumed no loss in precision 
from using the same value of b for both the A^^ and Ag 
distributions. The present procedure, however, clearly 
showed a small difference in the two, the b for the 
resonance distribution being less in every case, as it 
should, The work of Rush with Ra-Be neutrons [l6] showed 
this same tendency. Figure ^ is a plot of A^^ vs. r and 
Ag vs . r for the pure water run which includes the data 
computed from Equation (35) > normalized to the experimental 
data. Figure ^ also shows the marked deviation between 



3 ^ 



Activity (cpm) 




Figure ^ 
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the computed and experimental data for points near the 
source. This is to be expected since Equation (35) applies 
only at larger distances from the source. Similar plots 
for the poisoned water runs are presented in Appendix B. 

The values of b determined for all experimental runs are 
given in Table 3. 



TABLE 3 

RELAXATION LENGTH 



Boric acid 
(g/£) 


b(± 0.03 cm) 
Indium Resonance 
Activation 


b(± 0.02 cm) 
Total Activation 


_ _ _ 


10.30 


10.50 


1.16 ± .17 


10.25 


10 . il 5 


8.96 ±-.36 


10.25 


10.38 


27.81 ± .79 


10.22 


10.36 


^6.65 ± 1.2 


10.19 


10.35 


70.09 ± 1.3 


10.15 


10.35 


80.09 ± 1.9 


10.15 


10.35 



C. DETERMINATION OF AGE AND MIGRATION AREA 

A smoothed curve was drawn through data points less than 
r = 12.5 cm. on the A^^r and A^r versus r plots. For 
r > 12.5 cm. the data were extrapolated to infinity by means 
of a straight line having a slope equal to the negative 
reciprocal of the parameter b. The areas under the smoothed 
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curves were found by integrating the straight-line portions 

analytically and by using Simpson’s Rule for the remainder. 

Points from these curves were used to obtain the curves 

Acd^ Agr versus r. These latter curves vjere found 

to approach straight lines beyond r = ^5 cm. and areas were 

calculated accordingly. After the area under each of the 

_ 2 

four curves had been determined, the second moments rg and 

r^ were obtained from Equations (25) and (33). The neutron 

age to indium resonance and total migration area then 

followed from Equations (31) and (3^). The computer program 

used in these computations is Included in Appendix G. 

2 2 ^ H 

Figures 5 and 6 are plots of A^^r , Agr , A^j^r , and Agr 

2 2 

versus r for pure water. Plots of A^j^r and Agr versus 

r for the poisoned experiments are given in Appendix C. 

2 

Table h gives the values of and M,p determined for 

all experiments. The errors cited have been estimated from 
2 2 

Agjjr and A^r versus r plots by using the individual, 

unaveraged foil points to construct maximum and minimum 

2 

values for and Mrp [11] . 

A rough check on the experimentally determined values 

of age was made by following the procedure discussed in 

Section II-D with Equation (32). Values of activity as a 

function of r had been previously computed from Equation 

(35) in determining the closest value of relaxation length 

to fit the experimental points. This data was plotted on 

2 

semilog paper against r . A narrow band betv/een r = l8 
and r = 22 cm. presented the closest straight line 
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Figure 6 
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TABLE 



MIGRATION AREA AND NEUTRON AGE 



Boric Acid 
g/t 


M^^(cm^) 




T f 2 ^ 

Ij^Ccm ) 

(From Eq. ( 32 ) ) 


_ _ _ 


61.76 ± 


1 . 8^1 


52.19 ± 0.79 


51.75 


1.16 ± .17 


59 . ^1 ± 


0.68 


52.06 ± 2.3^ 


51.62 


8.96 ± . 36 


55.63 ± 


1.17 


52.02 ± 3.18 


51.62 


27.81 ± .79 


5 ^. 21 ± 


2.89 


51.96 ± 3.06 


51.55 


^6.65 ± 1.2 


53.^7 ± 


3.13 


51.29 ± 3.^2 


51.^7 


70.09 ±1.3 


53.37 ± 


1.57 


51.26 ± 3.11 


51.37 


80.09 ±1.9 


53.13 ± 


3.13 


51.20 ± ^.21 


51.37 



approximation. The slope of this line, equal to - 

^■^IN 

•from Equation ( 32 )^ was then determined by the method of 
least squares. A plot of A^^ vs. r^ for pure water with 
b = 10.3 is provided in Appendix D. Table Includes values 
of obtained in this manner for all experimental runs. 

The least square machine program used in this determination 
of age is Included in Appendix H. 

D. DETERMINATION OP DIFFUSION LENGTH 

Prior to determining Lrp from Equation ( 2 ^), It was 
necessary to apply a correction to to account for the 

additional age required to reach the thermal energy of 

0.025 eV. 






mm 



Estimates of this correction for hydrogeneous media are 

2 2 

found in the literature ranging form 0.3 cm [17] to 1.8 cm 
[11] , the latter value being that used by Valenti for 
measurements with a Pu-Be neutron source. In this investi- 
gation it was found that adding a correction of 1.8 cm to 
the value of determined for pure water yielded a value 

of 2.786 for Lrp. This is in excellent agreement with 
published results obtained by the more conventional method 
of measuring the exponential attenuation of thermal neutrons. 
It was also noted that 2.786 cm was in precise agreement 
with the value of L,p obtained from the empirical equation 
[ 1 ^] 



= 2.6^ + 0.0061 T cm 



(37) 



HoO 



for a temperature, T, of 2H°C. Accordingly, a correction 
of 1.8 cm for the age from indium resonance to thermal ener- 
gies was applied to the experimentally determined value of 
Tij^ for all runs and was then computed from Equation (24). 
Results are presented in Table 5* Errors have been estimated 
from a combination of those previously cited in Table 4 for 

^IN ^T * 



E. DETERMINATION OF THE THERMAL FLUX DISTRIBUTION ACCORDING 
TO TWO-GROUP THEORY 

A final check on the results obtained for pure H 2 O and 
the poison mixtures was made by computing the thermal flux 
distribution from Equation (15). The experimentally deter- 
mined values of L^ and Trp were used in the equation and the 
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TABLE 5 



MEASURED DIFFUSION LENGTH AT 24 °C IN 
PURE WATER AND 40°C IN POISONED V/ATER 



Boric Acid L^(cin) 

g/£ ' i 



1.I6 ± ,17 
8.96 ± .36 
27.81 ± .79 
46.65 ± 1*2 

70.09 ± 1.3 

80.09 ± 1.9 



2.7862 ± .14 
2.3559 ± .15 
1.3307 ± .14 
0.6664 ± .19 
0.6213 ± .20 
0.5594 ± .18 
0.3640 ± .30 



total macroscopic absorption cross-section of the poisoned 
medium was computed from 

r = z + 

“ “hydrogen Boron 

where Z was taken to be 

^hydrogen 

2 = Ro X 3a 

^hydrogen ^H 20 hydrogen 

and 

Z^ = N„ X a 

^Boron H^BO^ ^Boron 

The parameter a is the microscopic absorption cross-section 

^ 3 

and N„ p>o number of molecules of H^BO^ per cm 

li ^ 

determined from 



I 



'i 



N = 2 2 

Molecular weight 

3 

where p is the density in grams/cm and Aq is known as 

2 *5 

Avogadro ’ s number, which is 6.02^ x 10 molecules per gram 
mole. 

The computed values of flux were then normalized to the 
thermal saturated activity which was obtained by subtracting 
the covered foil from the bare foil data. The results were 
plotted together on semilog paper and compared. Figure 7 
gives this comparison plot for pure water. Similar plots 
for other runs are contained in Appendix E. Table 6 lists 
the values of determined for each run together with 
the maximum deviations found between the experimental data 
and flux distributions obtained from Equation (15)* 



TABLE 6 

COMPARISON OF EXPERIMENTAL RESULTS WITH SOLUTIONS 
0F THE TWO-GROUP THERMAL FLUX EQUATION 



Boric acid 
g/l 


f^(cm"^) 


MAX deviation 
{%) 


^ _ 


0.0222 


11.7 


1.16 ± .17 


0.0307 


lH.2 


8.96 ± . 36 


0.0882 


12.5 


27.81 ± .79 


0.2270 


7.2 


il6.65 ± 1.2 


0.3657 


3.5 


70.09 ±1.3 


0.5381] 


2.3 


80.09 ±1.9 


0.6120 


- - - 



^3 



I 




Figure 7 
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V. COMPARISON OF RESULTS WITH THEORY AND 
PREVIOUS EXPERIMENTAL MEASUREMENTS 



A. AGE AND THERMAL MIGRATION AREA 

Theoretical investigation of the age of fission neutrons 
in pure H2O has been extensively pursued and, for this case 
alone, there exist calculations of the actual slowing down 
distribution which can be compared in detail with measured 
distributions. Comparable results have not been obtained 
for non-fission sources due, primarily, to uncertainties in 
the energy spectrums. There is no account in the literature 
of an age calculation for a Pu-Be source. Coveyou and 
Sullivan [ 17 ] , however, have calculated the age to indium 
resonance, for monoenergetic sources in H2O using a 

Monte Carlo' age code which provides a rough comparison for 
this work. The plotted results of their calculations demon- 
strate a rapid rise in age with increasing source energy; 
and for an initial neutron energy corresponding to the 

average energy of a Pu-Be source ( 4.2 MeV) , these results 

2 2 

give a value for Tjj^ of 52 cm (t.5 cm ) . 

2 

Although clearly defined, thermal migration area, M^p , 

is not amenable to calculation because of theoretical 

difficulties encountered in describing the scattering of 

thermal neutrons [I8] . As was borne out by this experiment, 

however, thermal diffusion length is small in highly 

absorbing systems so that the predominant contribution to 
2 

M^ comes from the slowing down process. This allows some 
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( 



2 

comparison of measurements of and calculations of 

by treating the thermal diffusion as an additive correction. 

2 

The only measurements of and x^^^ from a Pu-Be source 

which have been published to date were those by Valenti 
and Sullivan for x^j^ [ 19 ], and by Valenti for M^^[ll]. 

Their results are compared with those obtained In this work 
in Table 7 . A factor which accounts for some of the dif- 
ference in these two measurements is a source effect cor- 
rection used by Valenti and Sullivan [ 19 ] . This Involved 
a correction to the indium resonance activity to account 
for the finite size of the source and was Intended to 
approximate a distribution closer to that expected from 
a point source. Due to the relatively low counting rates 
encountered in this work, it is estimated that a source 
effect correction would have Increased the activity of the 
foil nearest the source by less than \% [ 11 ] . 



TABLE 7 

COMPARISON OF RESULTS WITH PREVIOUS EXPERIMENTS 
Plutonium-Beryllium Neutrons in Pure Water 





Present Results 


Valenti & 
Sullivan [ 19 ] 


Valenti [ 11 ] 


Mrp^(cm^) 


61.76+ 1. 




61.07 


( 2v 

Tij^(cm ) 


52.19 ± 0.79 


52 . 8 ± 2.0 








and 



2 

In Figure 8, measured values of the parameters Mrp 

^IN plotted as a function of poison concentration. The 

differences between the two curves Is seen to reflect the 

1/v behavior of the boron absorption cross-section. At low 

2 

thermal energies absorption Is heavy and affects only M,p ; 

at energies above thermal, absorption falls off rapidly and 

2 

has a very small affect on both and 

B. THERMAL DIFFUSION LENGTH 

The close agreement between the theoretical value of Lrp 
In pure H 2 O at room temperature with that obtained In this 
experiment was pointed out In Section IV-D. 

Previously reported Investigations of Lrp dependence on 
the poison concentration of H^O were conducted using a ther- 
mal reactor as a neutron source [3, and a Sb-Be (25keV) 
source [5] . As a result there existed only a small com- 
ponent of fast-flux In the region of measurement which 
could easily be eliminated from the experimental data. 

This permitted relatively straightforward measurements of 
the exponential attenuation of thermal neutrons as discussed 
in Section II-C and, consequently, gave final results of 
high accuracy. 

Figure 9 compares the results obtained by Goddard and 
Johnson [3] with values of Lrp found in this experiment. The 
two sets of data are found to be In good agreement when the 
25°C temperatui'‘e difference is taken Into account. Since 
their results were obtained at 65°C, their values of Lrp 
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Poison Concentration (g/1) 




i\9 



Poison Concentration (g/1) 



should be higher according to Equation (37)- It can be seen 
from the curve representing the present measurements that 
absorbing solutions with concentrations of 70.09 and 8O.O9 
g/1 lay outside the Corngold limit (0.6 cm). 

C. TWO-GROUP THERMAL FLUX 

The thermal flux computed from Equation (15) is compared 
with measured values of saturated activity for each experi- 
mental run in Figure 7 and Appendix E. Due to the compli- 
cated manner in which Pu-Be neutrons slow down, these 
computations provided the only method of examining initial 
activity data for agreement with a distribution predicted 
by theory. It is evident from the maximum deviations given 
in Table 6 that the measured and calculated distributions 
come closer- in agreement with increasing poison concentration. 

The purpose of the experimental set up was to approxi- 
mate as closely as possible the "infinite medium" condition 
which was a basic assumption in the theoretical treatment 
of both age and diffusion length. That there was, in fact, 
no leakage in the axial direction is deemed a valid conclu- 
sion from the absence of detectible bare foil activity 
beyond 32 cm. from the source. An indication of the extent 
of lateral leakage was obtained by measuring total activity 
radially outward from the source in pure H 2 O. The extrapo- 
lated distance beyond the physical dimensions at which the 
flux can be assumed to fall to zero was found to be 1.29 cm. 
This result is displayed in Figure 10 from which it can be 
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seen that even in pure water the "infinite medium" assump- 
tion was not seriously violated. The minor differences 
between the measured and calculated flux distributions for 
concentrations of ^ 6.65 and 70.09 g/1 indicate the better- 
ment of the assumption at the heavier absorber concentration. 
This is most apparent from a study of the plots in Appendix 



VI. CONCLUSIONS 



The dependence of diffusion length on poison concentra- 
tion was found to agree well with that reported by Goddard 
and Johnson 13]. However, because diffusion length as 
obtained in this experiment is the square root of a small 
difference between two large numbers, the results given in 
Table 5 include sizable estimates of error. Errors in the 
measurem.ents from which diffusion length was obtained were 
due to limited experimental conditions, the most serious of 
which was the weak neutron source, resulting in a low activity 
and a small number of data points. This in turn Increases 
the statistical error. 

Previous measurements of neutron age and migration area 
from a Pu-Be source in pure water [11, 19] agree with the 
present results within the experimental errors. This report 
is the first of measurements of these parameters as a 
function of poison concentration. Accordingly, the results 
given in Table ^ should be of considerable value in extend- 
ing our understanding of the slowing down process. The 
errors associated with each measurement have been conserv- 
atively estimated and are meant to serve as bounds on the 
error, rather than a probable error. 

It is evident from the experimental data plotted for the 
heaviest concentration (80 g/1) that it was impossible to 
make a reliable age measurement in this case. The value 
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given for this run is therefore merely a logical estimate 
based heavily on the previous measurements. It follov/s from 
this that the numerical value listed for diffusion length 
at this concentration Is also unreliable. However, since 
it can be assumed that each of these values approximates 
the result one v;ould expect to obtain had a measurement been 
possible, they were used in Equation (15) to calculate a 
thermal flux distribution. The result of this calculation 
is plotted in Figure 11 together with the experimental data 
through which no curve could be drawn. Figure 11 vividly 
shows the irresolute behavior of the thermal distribution 
at the 80 g/1 boron concentration. It is concluded from 
this, and from the failure of the age measurement at this 
concentration, that the "maximum absorption" limit was 
exceeded in this experiment. The cadmium measurements of 
Goddard and Johnson together with the present results 
indicate that the diffusion length corresponding to this 
limit is about 0.5 cm, slightly less than the theoretical 
value. The absorption cross-section of the supporting 
medium at the time the final measurements were attempted 
was 0.6120 cm~^. 



5 ^ 



( 

t 



Thermal Activity (cpm) 




Figure 11 



APPENDIX A 



METHOD OF CHEMICAL ANALYSIS FOR DETERMINING 
BORIC ACID CONCENTRATION 

A, PREPARATION OF STANDARD 

A l-^gram sample of H^BO^ (Mol wt 6l.8^), finely ground, 
was dissolved In 300 ml of pure H2O . To this solution In a 
beaker were added 10 drops of phenolphthaleln indicator and 
25 ml of neutral glycerol. From a buret, standard NAOH was 
added dropwise until a distinct reddlsh-plnk color appeared. 
This marked the true end point for boric acid and established 
the molarity (Z) of the standard in moles/liter. 

B, ANALYSIS OF MODERATOR CONCENTRATION 

To 25 ml of moderator solution were added indicator and 
NAOH as described above. The standard alkali was then 
added dropwise until the true end point was reached. 

C , EXAMPLE 

1.33 ml of NAOH produced end point In solution 
containing 0.162 moles H^BO^. 

Z = “ 12.15 moles/liter standard alkali 

0,3 ml of standard alkali produced end point In 
25 ml of sam.ple drawn from moderator 

(25 ml)(molarity of sample) = (.3 ml)(Z) 

molarity of sample = .li ^6 moles/liter 

(61. 8 ^) (. 1 ^ 6 ) = 9*03 g /1 H^BO^ concentration of 
moderator ^ ^ 
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APPENDIX B 



PLOTS OF INDIUM RESONANCE AND TOTAL ACTIVITY VERSUS 
DISTANCE FROM THE SOURCE FOR POISON RUNS 
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Figure 13 




-Figure 



59 




Figure 15 



60 





Figure 16 
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Figure 17 
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APPENDIX C 




Figure 18 
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Figure 19 
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Figure 20 
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Figure 21 
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, Figure 23 
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APPENDIX D 



PLOT OF CALCULATED INDIUM RESONANCE ACTIVITY 
VERSUS p2 FOR RELAXATION LENGTH OF 10.3 cm 




Figure 2h 



APPENDIX E 



PLOTS OF THERMAL FLUX VS r FOR POISON RUNS 




Figure 25 
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Figure 26 
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Figure 27 
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APPENDIX P 



COMPUTER PROGRAM 1 



A, PURPOSE 

1. To correct counting data and compute the saturated 

activity in counts per minute and Ag), 

2. To compute the approximate mean counting rate and, 
from this, determine the standard deviation based 
on a normal distribution. 

2 2 

3. To compute the products A^^r and A^r . 



B. PROGRAM NOTATION 



D 

T 

POS 
C, CC 

ASAT, BSAT 

AMEAN 

W 

SIGMA 

P 



,116 



Decay constant for IN’ 

Counter resolving time 
Distance from source 
Ten minute counts (side 1, side 2) 
Saturated activity (side 1, side 2) 
Approximate mean counting rate 
Correction for foil weight 
Standard deviation 
Cadmium correction factor 
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APPENDIX F 
COMPUTER PROGRAM 1 



IMPLICIT REAL*8( A-H,C-Z) 

WRITE (6,1) 

1 FORMAT! T5, '01 STANCE • ,T25 SAT ACTI VI TY • ,T^6 » • ME AN ' , 
1T87, 'SIGMA ' //) 

D=1 .280-2 
T=6.C0-6 
D050C M=l,2 
0050 M=l,9 

READ(5t2) POS,C ,X,Y,Z,CC,YY,ZZ 

2 FORMAT! 6F9.0) 

READ! 5, 6) W 

6 FORMAT! 1F9.0) 

DI ST = P0S’!‘0.25AD01 
C=W«C/! l.COO-C*T) 

CC = W^CC/! 1.0D0-CC*T) 

ASAT = D’i'C/! ! 1 . 00 C-DE XP ! -D=«=X ) ) ’'^ ! OF XP ! -D*Y ) - DFX P ! -0* Z ) ) ) 
BSAT=D>yCC/! ! 1 . 00 0-DE XP ! -D^^X ) ) -i' ! DEX P ! -OvYY ) -CcX P !-0* Z Z ) 

1 ) ) 

SAT=!ASAT+DSAT) /2.000 
AHEAN=!C+CC) /2.0001 
SIGMA=AKEAN*«0. 500 
F=1.1100 
FSA T=F=^SAT 

WRITE! 6,4 ) 01 ST, FSA T, A MEAN, SIGMA 

4 FORMAT! T6,F5.2,6X,025.16,6X,025.16,5X,D25.16) 

A=FSAT>f'DI ST*«2. ODO 

8=FSATvOI ST*«4. ODO 

8 FORMAT! / , T1 2 , F 5 . 2 , 6X ,02 5 . 1 6 ,6 X , 025 . 1 6 ) 

WRITE! 6,8) DIST,A,B 

50 CONTINLE 
0060 N=l,9 

R£AD!5,5) POS ,C , X , Y , Z , CC , YY , ZZ 

5 FORMAT! 8F9.0) 

RcAO!5,7) W 

7 FORMAT! 1F9.0) 

OIST=P0S*0. 254001 
C=W>^C/! 1.000- C’i'T) 

CC=W’i'CC/! 1 .0rj0-CC*T) 

ASAT=O^C/! ! 1. 0DC-DEXF!-D*X) ) * ! DEXP ! -O^Y ) - DEXP ! -D* Z ) ) ) 
BSAT = DvCC/! ! 1.00 0- 0EXP!-0>«'X) ) * ! DEX P ! -0«YY ) -DEX P ! -D* ZZ ) 

1 ) ) 

SAT=! ASAT+BSAT) /2. ODO 

AM,EAN=!C+CC)/2.0D01 

SIGMA=AMEAN>f--^0. 500 

WRITE! 6,3) DIST, SAT,AMFAN,SIGMA 

3 F0RMAT!T6,F5.2 ,6X,D25.16 , 6 X ,02 5 . 15 , 6 X , 025 . 1 5 ) 

A=SAT«DI ST*y2. ODO 

B = SAT*DIST’^--.'4.0D0 
WRITE! 6,9) 0IST,A,B 

9 FORMAT! / ,T12,F5.2 ,6X,D25.16 ,6X,D25. 16) 

60 CONTINUE 

500 CONTINLE 
STOP 
END 
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APPENDIX G 



COMPUTER PROGRAM 2 



A. PURPOSE 

1. To compute the area under curved portions of the 

2 2^ H 

curves A^^r , A^r , A^^r and A^r versus r using 
Simpson's Rule. 

2. To compute area under straight line portions of 
these curves analytically. 

2 

3. To compute total areas and, from this, M,p and 

using Equations (13) and (l8). 

To compute by correcting for additional age 

required to reach thermal energies. 

5. To compute diffusion length by solving 

Lrj. = (M^2 - 

B. PROGRAM NOTATION 

2 2 ^ 

y(I) Ordinate values of A^^r , A^r , A^^r and 

Agr*! 

H Abscissa Interval lengths in cm. used to 

obtain y(I) values. 

N Number of abscissa intervals, H. 

S(J) Simpson area approximations 

B(I) Relaxation Length values 

2 2 4 

Z(I) Ordinate values of A^^r , A^r , A^j^r 

and Agr^ corresponding to abscissa values 
of r at which distributions become exponen- 
tial . 
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i 

\ 

I 



i 




Ad) 

AA(I) 

MIG(I) 

AGE(I) 

Y(I) 



Area under exponential portion of curves. 
Total area under curves. AA( I) = S( J) + A(I) . 
Thermal Migration area values. 

Values of neutron age to Indium Resonance. 
Diffusion Lengths. 
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APPENDIX G 
COMPUTER PROGRAM 2 



COMPUTING SATURATED ACTIVITY TIMES DISTANCE SQUARED AREAS 
IMPLICIT REAL*8( A-H ,C-7) 

R£AL^8MIG 

DIMENSION Y( 80) , S( 28) ,B( 28) , 7 ( 28 ) , A ( 23 ) , AA ( 2 3 ) , 
1MIG( 28 ) ,AGE( 28) 

DO 80 J=1,1A 

READ! 5 ,10) ( Y( I ) ,I =1 ,11) 

10 FORMAT! IFIO.O) 

F=1.25 

N=10 

SUM2=0.CPC 
005C L=2,N,2 

50 SUM2= SLM2+Y(L) 

SUM3=0.CDC 
M=N- 1 

D060 K=2,M,2 
6C SUM3= SIM3+Y( K ) 

S( J )=H/0. 3D1*( Y( 1 ) +0.ADl’!‘SUM2 + 0. 2D1*SUM3+Y ( N + 1 ) ) 
WRITE(£,7C) S(J) 

7C FORMAT! Tl^, • SI MPSON A P PRCX I M AT I CN= ' , D2 5 . 16/ / / ) 

8C CONTINUE 

COMPUTING ACTIVITY TIMES DISTANCE TO THE FOURTH ARFAS 
D090 J=15,28 

READ!5,11) ! Y! I ) ,I =1 ,19) 

11 FORMAT! IFIO.O) 

H=2.5DC 

N=18 

SUM2=0.0D0 
D051 L=2,N,2 

51 SUM2=SLM2+Y! L) 

SUM3=0.CD0 

M=N-1 

D061 K=3,M,2 
61 SUM3=SLM3 + Y! K) 

S(J)=H/C.3D1»!Y!1) +0 . 4D1 >?^SUM2+0. 201>i'SUM3+Y ( N+1 ) ) 
WRITE!6,71) S!J) 

71 format (T i 9, • SI MP SON' ,025.16///) 

90 CONTINUE 

COMPUTING AREAS UNDER LINEAR PORTION OF CURVES 
READ!5,I) !B! I ) ,Z! I ) ,1=1 ,23) 

1 FORMAT! 2F1C. 0) 

DO 2 1 = 1,26 

A! I ) = Z! I )/B! I ) 

WRITE! 6,3) R! I ) , A! I ) 

3 FORMAT! 2E1A. 7///) 

2 CONTINUE 
00100 1=1,28 

100 AA! I ) = A! I ) + Sl I ) 

D0200 1=1,13,2 

MIG! I )=AA{ I + IA) /! AA!I 1^6. 0) 

WRITE! 6,201) MIG! I ) 

201 FORMAT! T19 , 'MIGRATION ARE A= • , 02 5 . 1 6 /// ) 

2CC CONTINUE 

D0300 1=2, U, 2 

AGE! I )=AA! I + 1A)/!AA!I )=!'6.C) 

WRITE! 6,301) AGE! I ) 

3C1 F0RMAT1T1?,'AG£=',D25.16///) 

20G CONTINUE 

D0600 1=2, lA, 2 

Y!I)=!MIGIJ-1)-!AG5!I)+1.3))**0.5 
WRITF!6,AC1) Y! I ) 

ACl FORMAT! T19, 'DIFFUSION LE KG TH= ' , D25 . I 6/ / ) 

ACC CONTINUE 
STOP 
END 
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APPENDIX H 



COMPUTER PROGRAM 3 



A. PURPOSE 

2 

1. To determine slope of curve ^nA^^ vs. r In the 
range l8.iJ £ r £ 22 cm. by the method of least 
squares . 

2. To compute from the relation 

, 1 
slope = - ]T- — 



B. PROGRAM NOTATION 

F(I) Relaxation length values 

G(I) Values of neutron age and thermal migration 

area 



X(I) 

Z(J) 



Y(I) 

SIMQ 



a 



2 



Distances from the source 
Solutions to equation 




Natural logarithm of Z values 
Y = £n Z = a^ + a 2 X 

Subroutine for solving system of equations 
CA = B: 





C 


A 


= 


B 




7 








7 


7 


2 X. 

i=l ^ 




^1 




^ ^i 
1 = 1 ^ 


7 


7 2 

X X X 

1=1 _ 








7 


5: 

i=l 




^2 




I X y 
1=1 ^ ^ 



Slope of straight line solution equal 
to the negative reciprocal of 
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APPENDIX II 
COMPUTER PROGRAM 3 



DIMENSION X( 9) ,Y{9) ,7(9) ,XC(2) ,C(2,2) , 

1B(2) ,A( 2) ,D( 7) ,F (7) ,G(9) 

RcAO(5,5) (D( I ) ,F( I ) ,I =i ,6) 

5 FORMAT! 2F9.0) 

READ(5,5t) (G( I) ,1=1,6) 

56 F0RMAT(F9.0) 

D05 L=l,6 
X( 1) = 1E.A1 
X( 2)=ie.7 
X( 3)=19.2 
X ( A ) = 1 9 . 5 
X( 5)=19.7 
X( 6)=2C.O 
X{ 7)=22.C 
DD9 J=l,7 

Z( J ) = ( D( L )«eXP( - X( J ) /F (L) ) ) /X( J) *^>^2 .0 
9 CONTINUE 
Don 1=1,7 

11 Y( I )=ALOG(Z( I) ) 

DO 12 1=1,7 

12 X( I ) = X( I )««2.0 
SUMMING POWERS OF X( I ) 

00 20 K=l,2 
XC(K )=C.O 
CO 15 1 = 1,7 
15 XC(K) = XC(K)+X( I 
20 CONTINUE 
SUMM ING Y( I ) 

YC=C.O 
CO 22 1=1,7 
22 YC = YC + Y( I ) 

YX=0.0 
D025 1=1,7 
25 YX=YX+ Y( I )*X( I ) 

FORMING TWO EQUATIONS CA=0 WITH ELEMENTS IN COLUMNS 
DOAO 1=1,2 
D035 J=l,2 
JUGl=I+J-2 
IF( JUGl )33,33,31 
33 C( 1, 1)=7.C 
GOTO 35 

31 C( I, J)=XC( JUGl) 

35 CONTINUE 
40 CONTINUE 
B( 1)=YC 
B( 2) = YX 

CALL SIMQ(C,B,2,KS) 

WRITE(6,5C) KS 
50 FORMAT! T5, 'ERROR = ' ,13) 

DO 52 1 = 1,2 
52 A! I )=R! I ) 

WRITE!6,54) A!l) ,A!2) 

54 FORMAT! T5, 'A!1)=',E14.7/T5,'A!2)=',E14.7) 
AGE=0. 2S/A! 2) 

WRITE!6,55) AGE 

55 FORMAT! Ti5, *AG5= ' ,E14. 7//) 

6 CONTINUE 
STOP 

END 
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APPENDIX I 



RESULTS OF COMPUTER PROGRAM 1 
TABLE 8 

PURE WATER; INDIUM RESONANCE ACTIVATION 



r(cni) 


Acd(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


298.0691 


233.2133 


15.2713 


9.52 


20 i|.l 689 


157.3879 


12.5^5^ 


10.79 


155.0969 


118.2727 


10.8753 


12.07 


12^1.0281 


89.9505 


9.^8^^2 


13 . 3 ^ 


78 .ill 97 


19.3193 




1^.61 


67.7517 


50. il2il5 


7.1010 


15.88 


50.9003 


39.3977 


6.2767 


17 . 1 ^ 


H 2 . 2 l\ 0 H 


32.7216 


5.7203 


18. ill 


3^<.2155 


25.6226 


5.0618 



TABLE 9 

PURE WATER; TOTAL ACTIVATION 


r ( cm) 


AgCCPM) 


Mean Count 


Standard 

Deviation 


8.26 


i|il 59.2573 


3663.2972 


60.5252 


C\J 


35 il 9.0267 


29l8.2ill8 


5 i |.0207 


10.79 


27 il 7 . 3121 


2259.8505 


il 7.5378 


12.07 


1930 . 577 il 


l66i| . 8222 


i|0. 8022 


13.3^ 


Iil 37 .il 071 


Il 8 l. 5 il 08 


3 ii .3735 


l^\. 6 l 


1090.0382 


900. Iill8 


30.0023 


15.88 


852 . 3 il 89 


692. 328i| 


26.3121 


17.1^ 


6 il 2.0987 


529.1137 


23 . 002 i| 


18. ill 


5 ill .2537 


i|il 9.0388 


21.1905 
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TABLE 10 



H^BO^ CONCENTRATION l.l 6 g/ 1 ; INDIUM RESONANCE ACTIVATION 



r (cm) 


Acp(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


503. 1124 


379.3378 


19.4766 


9.52 


322.7215 


240.9674 


15.5231 


10.79 


233.0711 


176.0013 


13.2665 


12.07 


233.6401 


175.2742 


13.2391 


13 . 3 ^ 


172.9660 


129.3376 


11.3726 


li^.61 


116.5811 


90.3156 


9.5035 


15.88 


74.4533 


57.2474 


7.5662 


17 . 1 ^ 


94.6393 


72.5835 


8.5196 


l8.i<l 


48.3223 


36.5259 


6.0436 



TABLE 11 



H^BO^ CONCENTRATION l.l 6 g/lj TOTAL ACTIVATION 



r( cm) 


Ag(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


6261.3213 


5172.9621 


71.9233 


9.52 


4148.3711 


3446.4026 


58.7061 


10.79 


3113.6634 


2596.8206 


50.9590 


12.07 


2229.4831 


1854.0849 


43.0591 


13.34 


1681.5223 


1400.5575 


37.4240 


14.61 


1265.8901 


1080.4300 


32.8698 


15.88 


986.9352 


842.4719 


29.0254 


17.14 


768.5511 


639.9204 


25.2966 


18.41 


635.0632 


527.7004 


22.9717 
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TABLE 12 



H^BO^ CONCENTRATION 8.96 g/ 1 ; INDIUM RESONANCE ACTIVATION 



r (cm) 


AcpCCPM) 


Mean Count 


Standard 

Deviation 


8.26 


il 09 . 9^67 


31^.5537 


17.7356 


9.52 


291.^^77 


213.2623 


lil.6035 


10.79 


229.5099 


167.0131 


12.9233 


12.07 


167.2182 


123.1536 


11.0975 


13.3^ 


lil 2.2078 


103.2886 


10.1631 


1^.61 


105.7302 


81.3866 


9.0215 


15.88 


67.2199 


51 .ill 56 


7.1705 


17 . 1^1 


83.3672 


62.9565 


7.93^5 


18 . ill 


ilil . il 031 


33 . 52 il 5 


5.7900 



TABLE 13 

H^BO^ CONCENTRATION ^8.96 g/ 1 ; TOTAL ACTIVATION 



r (cm) 


Ag(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


2713.8679 


2326.7305 


48.2361 


9.52 


1704.0191 


1413.4715 


37.5961 


10.79 


1188.5005 


986.2288 


31.4042 


12.07 


871.8115 


725.4009 


26.9333 


13 . 3 ^ 


639.1190 


530.0796 


23.0235 


lil. 6 l 


487.743 


405.6754 


20 . 14 i 4 


15.88 


373.6443 


310.4774 


17.6204 


17.1^ 


279.7307 


226.8509 


15.0616 


18 . ill 


269.3002 


220.0322 


14.8334 



w 





TABLE Ik 



H^BO^ CONCENTRATION 27.81 g/ 1 ; INDIUM RESONANCE ACTIVATION 



r ( cm) 


A^p(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


285.27^7 


218.7960 


liJ.7918 


9.52 


236.3981 


183 . 69^^9 


13.553^ 


10.79 


ns.okik 


138.1119 


11.7521 


12.07 


135.2377 


103.^736 


10.1722 


13 . 3^1 


lii 6 . 917 ^ 


107.5823 


10.3722 


1^.61 


89.3376 


66.5^02 


8.5721 


15.88 


69.33^3 


51.0123 


7.1^23 


17,1^ 


93.919^ 


73.3271 


8.5631 


18 . ill 


i< 2.0807 


32.5^63 


5.7049 



TABLE 15 



H3BO3 


CONCENTRATION 


27.81 g/ 1 ; TOTAL 


ACTIVATION 


r (cm) 


Ag(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


1313.3395 


1118.5688 


33.4450 


9.52 


855.6882 


729.7254 


27.0134 


10.79 


658.1651 


558.1721 


23.6256 


12.07 


503.8277 


427.9465 


20.6868 


13.34 


360.5809 


306.2933 


17.5012 


14.61 


399.0305 


332.2687 


18.2282 


15.88 


206.7903 


176.8829 


13.2997 


17.14 


186.5858 


158.2384 


12.5793 


18.41 


129.9547 


110.3157 


10.5031 
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TABLE 16 



H^BO^ CONCENTRATION ^ 6.65 g/ 1 ; INDIUM RESONANCE ACTIVATION 



r ( cm) 


A^p(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


383 . 1^35 


297.7103 


17.2543 


9.52 


282.^1397 


212 . 2^161 


14.5761 


10.79 


207.^628 


153.580^1 


12.3927 


12.07 


162.067^ 


119.95^9 


10.9524 




150.3377 


113.82^15 


10.6688 


1 ^ 1.61 


90.8321 


67.9496 


8.2432 


15.88 


95.1983 


71.8575 


8.4768 


17 . 1 ^ 


5^.^1682 


40.6000 


6.3718 


18 . ill 


79.3777 


57.6782 


7.5946 



TABLE 17 



H^BO^ CONCENTRATION ^l 6.65 g/ 1 ; TOTAL ACTIVATION 



r ( cm) 


A3(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


1006.1537 


835.3732 


28.9028 


9.52 


652.8737 


542.5276 


23.2922 


10.79 


455.5653 


377.8619 


19.4387 


12.07 


364.3542 


304.2721 


17.4434 


13.34 


310.0154 


251.2099 


15.8496 


14.61 


229.8385 


186.2102 


13.6458 


15.88 


199.8872 


159.1900 


12.6172 


17.14 


110.6291 


91.8081 


9.5816 


18.41 


88.0501 


73.0562 


8.5473 
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TABLE 18 



H^BO^ CONCENTRATION 70.09 g/ 1 ; INDIUM RESONANCE ACTIVATION 



r ( cm) 


A^p(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


379.0630 


285.6253 


16.9004 


9.52 


282.1943 


214.0142 


14.6292 


10.79 


230. 4439 


165.9322 


12.8814 


12.07 


160.5295 


120.9588 


10.9981 


13.34 


106.3356 


79.5638 


8.9198 


14.61 


95.0191 


69.8891 


8.3599 


15.88 


72 . 4740 


55.6762 


7.4617 


17.14 


51.1528 


38.2288 


6.1829 


18.41 


54.6934 


41.7976 


6.4651 



TABLE 19 

H^BO^ CONCENTRATION 70.09 g/ 1 ; TOTAL ACTIVATION 



r(cm) 


Ag(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


824.4368 


683.2489 


26.1391 


9.52 


558.2886 


462.7575 


21.5117 


10.79 


426.2032 


354.2170 


18.8206 


12.07 


316.5702 


261.6614 


16.1759 


13.34 


228.9763 


189.6963 


13.7730 


14.61 


192.6427 


159.1687 


12.6162 


15.88 


137.7566 


114.1718 


10.6851 


17.14 


93.6648 


79.3185 


8.9061 


18.41 


72.2664 


60.5687 


7.7826 
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TABLE 20 



H^BO^ CONCENTRATION 80.09 g/ 1 ; INDIUM RESONANCE ACTIVATION 



r ( cm) 


Acd(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


ill4. 8513 


307.2195 


17.5277 


9.52 


256.65^1 


196.2806 


14.0100 


10.79 


239.3121 


177.5645 


13.3253 


12.07 


I 6 l. 0 i <97 


122.6576 


11.0751 


13 . 3 ^ 


113 . 0 i <85 


87.4346 


9.3506 


li|.6l 


99. 86i<7 


77.9423 


8.8285 


15.88 


90.3185 


68.9324 


8.3026 


17 . 1 ^ 


SQ . kk^l 


45.0622 


6.7128 


18 . ill 


35.^877 


27.8769 


5.2799 



TABLE 21 



H3BO3 


CONCENTRATION 


80.09 g/1; TOTAL 


ACTIVATION 


r (cm) 


Ag(CPM) 


Mean Count 


Standard 

Deviation 


8.26 


743.8847 


614.4305 


24.7877 


9.52 


494.7427 


419.9927 


20.4937 


10.79 


360.5303 


309.5135 


17.5929 


12.07 


302.0725 


261.4078 


16.1681 


13.34 


212.3189 


175.3501 


13.2419 


1 — 1 
VT) 

rH 


149.1620 


128.0910 


11.3177 


15.88 


127.8776 


105.3823 


10.2656 


17.14 


116.5993 


99.1096 


9.9554 


18.41 


70.7709 


53.4546 


7.3113 
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was exceeded in this experiment. The absorption cross- 
section of the supporting medium corresponding to the 80.09 
g/1 concentration was 0.6120 cm~l. The present results,^ 
together with previous experimental data, provide an esti- 
mate of diffusion length at the absorption limit of 0.5 cm, 
Sj.ightly less than that corresponding to the theoretical 
Corngold limit (0.6 cm). 
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